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> I ■ Charge balance functions provide insight into critical issues concerning hadronization and trans- 

it ' port in heavy-ion collisions by statistically isolating charge/anti-charge pairs which are correlated by 

charge conservation. However, distortions from residual interactions and unbalanced charges cloud 
the observable. Within the context of simple models, the significance of these effects is studied by 
constructing balance functions in both relative rapidity and invariant relative momentum. Methods 
are presented for eliminating or accounting for these distortions. 

T— I ■ 
> 

in 

^ , I. INTRODUCTION 

o 

■ Charge balance functions were suggested as a means for addressing fundamental questions concerning hadronization 
' in relativistic heavy ion collisions The most pressing issue concerns whether hadronization is delayed in such 
, reactions beyond the characteristic time scale of 1 fm/c, i.e., is a new phase of matter created? A delayed hadronization 
of a gluon-rich medium would mean that many charge- ant icharge pairs would be created late in the reaction and then 
be more tightly correlated to one another in momentum space. Charge balance functions are designed to identify such 
' I ■ charge /ant icharge pairs on a statistical basis. Unfortunately, the ability to identify balancing partners is compromised 
O [ by two effects. First, surplus charge, originating from the non-zero baryon number and charge of the colliding nuclei, 
^ ■ pollutes the balance function. Secondly, interactions of a balancing pair with the other charges effectively polarize 
the other particles and distort the shape of the balance function. In this paper, the magnitude of such distortions are 
exhibited within the context of simple blast- wave models, and means for eliminating or reducing these distortions are 
presented. 

Charge balance functions are based on conditional distributions. 
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o.olo^ - 1 f iV+-(Pi,P2)-iV++(Pi,P2) . iV-+(Pl,P2)~A^— (Pi,P2) l ... 
^^^'^^'> = 2 1 N^) + Nli^) / ■ 

Here, Nab{Pi, P2) counts the incidences where a particle of charge a is observed with momentum in a region defined 
by Pi while a particle of charge h is observed that satisfies the momentum constraint P2 . The second constraint could 
be any function of the momenta of the two particles, e.g., the relative rapidity. Put into words, the balance function 
measures the chance of observing an extra particle of opposite charge given the observation of the first particle. Balance 
functions are related to charge fluctuations which can be used to investigate similar issues [3. IsL 13. IsL la. IsL 1^ llollTl| . 
The advantage of balance functions is that they represent a more differential measure. 

For a neutral system, every charge has an opposite balancing charge and the balance function would integrate to 
unity. 

Y,B{P2\Pi)^l. (2) 

The normalization is reduced if not all particles carrying the charge are included, e.g. only 7r+ and n~ are evaluated, 
thus neglecting the chance that the electric charge is balanced by a kaon or a baryon, or that the detector has less than 
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a perfect acceptance. If P2 refers to the relative rapidity, and Pi includes all measured particles, B{P2 — AF) provides 
the probability that a balancing charge was observed with relative rapidity ISY . Since much of the charge observed 
in a relativistic heavy ion collision should be produced at hadronization, a delayed hadronization should result in a 
tighter space-time correlation between balancing charges. Due to the large collective flow fields in these reactions, a 
tighter correlation in space-time translates into a tighter correlation between the final momenta. Therefore, a delayed 
hadronization should be signaled by a narrower balance function when plotted as a function of relative momentum or 
relative rapidity. 

One of the most enticing results from RHIC is the preliminary measurement of balance functions by the STAR 
collaboration In accordance with expectations for delayed hadronization, the balance functions appear to nar- 
row with increasing centrality of the collision. However, given the nascent stage of these observations and of the 
phenomenology, it should be emphasized that numerous questions remain concerning the interpretation of such a 
measurement. To that end, several issues were pursued in a previous paper, including the effects of Hanbury-Brown 
Twiss correlations, detector acceptance and the relation to charge fluctuations 

In the same spirit as that paper, more issues will be addressed in this study. In the next section, the benefits 
analyzing balance functions in other observables, e.g, the invariant relative momentum, will be addressed. In addition 
to allowing one to analyze the contribution from specific resonances, it will be shown that such observables help clarify 
other issues such as the interplay of collective flow and cooling. 

Balance function analyses are based on the assumption that all charges have balancing partners. This is not true in 
relativistic heavy ion collisions due to the presence of the initial protons and neutrons which bring about an imbalance 
of baryon number, electric charge and isospin. In section UlTl the distorting influence of the surplus positive charge is 
investigated and a modified balance function observable is proposed that would eliminate such effects. 

The subsequent section contains a detailed study of the effects of inter-pair correlations. By extending the model 
presented in [3| to balance functions in Qinv, it appears that the Hanbury-Brown Twiss (HBT) correlations cause a 
more noticeable distortion, especially in the most central collisions. The source of these residual effects is analyzed 
in detail, and the degree to which these distortions can be accounted for is discussed. The final section presents a 
summary of what further work must be done in analyzing and interpreting this class of observables. 



II. ANALYZING THE BALANCE FUNCTION IN Qinv 

In reference [l| balance functions were evaluated as a function of relative rapidity. Like two-particle correlation 
functions, the balance function is a six-dimensional quantity and new insights can be gained by performing different 
cuts or binnings. Specifically, we focus on performing analyses in terms of the invariant relative momentum, i.e. the 
relative momentum as measured by an observer moving with the velocity of the two-particle center of mass. We find 
that these variables yield clearer insight for interpreting the physics of the balancing charges, as well as providing a 
better illumination of the distorting effects which are the subject of this study. 

The relative momentum of the two particles is defined, 

qa = [Pa,a - Pb,a) - Pa (3) 

{.Pa.a Ph-Oi) -Pa • 



Here, the total momentum of the pair is P and the center-of-mass energy of the pair is ^/s = (pa + Pb)'^- For two 
particles of the same mass, the last term can be neglected. The invariant momentum is then 

Qinv-^q --[Pa-Pb) H ^2 ■ W 

For pion-correlation studies, it is conventional to define three projections of the relative momentum, Qiong, Qout and 
Qside [l^lT^ll5l |. These components measure the projections of q along the beam axis, the outwards direction (defined 
by the pair's transverse momentum) and the sidewards direction (perpendicular to the pair's transverse momentum 
and to the beam axis). Motivated by the semi-boost invariant nature of the collision geometry at RHIC, Qiong is 
usually measured in a reference frame moving with the beam velocity of the pair. Although not typically invoked in 
correlations studies, one can also perform a second outwards boost to a frame where the total spatial momentum of 
the pair is zero. In this frame the three components, Qout, Qsidc and Qiong sum to Qinv 

Qlong + Qsidc + Qout — Qinv (5) 



In terms of laboratory momenta P and q, these components are: 
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Qlong , „ {PoQz - Pzqo) (6) 
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Here, Pt ~ [{Pa,x +Pb,x)'^ + {Pa,y ''rPb.y)]^^^ IS the transverse momentum of the pah. These components differ from the 
common convention for HBT in that Qout is defined as the relative momentum in the pair frame, whereas in HBT the 
usual convention is to ignore the second boost which means that the three components do not sum to Qinv In fact, 
the factor {s + P^) / s in the definition of Qout is simply the Lorentz gamma factor corresponding to the transverse 
boost to the two-particle rest frame. 

Analyzing balance function in terms of Qmv simplifies interpretation with thermal models by eliminating the 
sensitivity to collective flow. Blast wave models are based on thermal emission from sources which move to account 
for the collective flow of the exploding matter. Collective flow affects the spectra, but leaves the invariant momentum 
differences unchanged if the two particles originate from the same space-time point of the blast wave. Hence, plotting 
the balance function in invariant momentum variables would minimize the confusion associated with the collective 
flow as the width would only depend on the local thermal properties of the individual sources. If a particle and its 
balancing particle were always emitted close to one another in coordinate space, the width of the balance function 
would principally be a function of the breakup temperature with no sensitivity to collective flow, assuming a uniform 
detector acceptance. 

To illustrate the complications of using rapidity differences rather than Qinv: one may consider a thermal source 
where the width in Qiong is determined by the temperature. The separation of two tracks in rapidity is then, 

Ay-^, (7) 

nit 

where mt is the transverse mass of the particles. Since collective flow affects the distribution of transverse masses, the 
balance function widths for localized thermal sources would depend on the collective flow in the data when plotted 
in relative rapidity. Although it is easy to account for collective flow in a theoretical model, the interpretation of 
experimental results is simplifled by performing the analysis in Qinv 

Furthermore, assuming thermal emission with highly localized charge conservation, the balance function would be 
isotropic with respect to the direction of the relative momentum, e.g., the width in Qside would equal the width in 
Qiong- For early production of charge, one expects string dynamics or diffusion to lead to an anisotropic balance 
function as the balancing charges should separate signiflcantly in coordinate space along the beam axis due to the 
extremely large velocity gradient along the beam axis at early times, dvz/dz = l/r. Thus, in addition to the width 
of the balance function in Qinv, the behavior of the anisotropy as a function of the collision's centrality provides a 
crucial test of the mechanism for charge creation and transport. 

To illustrate the sensitivity of a balance function in terms of these analyses, we consider a simple blast wave model 
where the collective transverse motion is assumed to rise linearly with the radius. Of the numerous parameterizations 
of the blast wave model, it is assumed that the sources have transverse rapidities governed by a simple distribution, 

dN _ I constant, yt < tanh"^(winax) ^-g-j 



ytdyt \ 0, yt> tanh ^{v 



'max } 

Here, yt is the transverse rapidity, tanh{yt) — vj_. The distribution of longitudinal rapidities is assumed to be uniform. 
For our calculation we assume that these sources emit isotropically in the source frame according to a temperature 
T. A balancing positive and negative pion are assumed to be emitted from sources with the same longitudinal 
and transverse rapidity. Figure ^ illustrates the sensitivity to the temperature by presenting balance functions for 
three temperatures, 90 MeV, 120 MeV and 150 MeV. The balance function is clearly narrower for lower breakup 
temperatures. The results of Fig. ^ are insensitive to the choice of Umax- However, it must be stressed that the 
sensitivity would return if the balance function is analyzed in a finite acceptance. 

The balance function could also be binned in any of the three projections, Qiong, Qout and Qsido, rather than in 
Qinv If the balancing pairs were to always originate from sources with the same collective velocity, the balance 
function would be identical in all three variables. However, if the balancing particles were to diffuse relative to one 
another, the shape of the balance function might become decidedly non- isotropic. For instance, if charge is created 
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early in a RHIC collision, the balancing charges might easily separate along the beam axis and ultimately be emitted 
from regions with different rapidities. Figure [5] presents the widths of balance functions assuming that the balancing 
particles independently dissipated and were each ultimately emitted with sources moving with a spread of rapidities 
characterized by cr^. A Gaussian form for the diffusion was assumed for the distribution of source rapidities, Us- 

P(y.)-exp(^). (9) 

This extends the distribution of Qiong while leaving the distribution of Qsidc unaffected and the distribution of Qout 
only slightly affected by boost effects. As can be seen in Fig. |3 the disparity in the three widths should be easily 
observed for this example where the temperature was chosen to be 120 MeV and the radial collective velocities were 
between zero and Wmax = 0.7c. 

It should be difficult to discern the difference between thermal broadening and dissipation of balancing charges into 
regions with different collective flow. However, other observables provide insight into breakup temperature, mainly the 
comparison of proton and pion spectra [l^ . Once one knows the breakup temperature, it is possible to fit parameters 
that describe the diffusive spread, e.g. ct^. Furthermore, a thermal fit to data where the diffusive terms are set to 
zero provides an upper bound for the breakup temperature. 

For the reasons above, much of the analysis of the following sections will be given in Qinv An additional advantage 
of using Qinv is that it allows one to identify the contributions from specific resonances which contribute peaks to the 
balance function when plotted in Qinv It is our hope that experimental analyses will also switch to these variables. 



III. THE EFFECTS OF SURPLUS POSITIVE CHARGE 



Not all charges have balancing partners. In a Au+Au collision at RHIC, the two gold nuclei provide 158 unbalanced 
protons and 236 unbalanced neutrons. These pollute the balance function by providing unbalanced electric charge, 
baryon number and isospin. For detectors like STAR, these effects are lessened by the fact that most of the surplus 
charge is at high rapidity and outside the experimental acceptance. However, the effect should become more significant 
if the balance function is constructed for a set of charges, e.g. pp, for which there is a significant imbalance of one 
charge vs. the opposite charge. Our goal in this section is to offer a revised procedure for producing balance functions 
from data that would subtract the pollution due to the surplus charge. More precisely, we wish to define a balance 
function that would ignore any additional unbalanced charges that are not correlated with one another or with 
pair-wise created charges. 

In order to demonstrate the effects of the polluting surplus charge, we introduce a notation where distributions, N, 
count charges which are divided into three categories. The subscripts "+" and "-" will refer to positive and negative 
charges which are created in tandem. The subscripts "(5" will denote the surplus positive charge. The balance function 
will be re-evaluated after inserting the replacements, 

N+-^N++Ns, 

- - - - (10) 

N++-^N++ + Ns+ + N+s + Nss, 

into Eq. iQJ. 

^ 'I 2iN+iP^) + NsiPi)) N+iPi) 

Ns{Pi)N+s{Pi,P2) - N+{Pi)Nss{Pi,P2) 
2{N+{Pi) + Ns{Pi))N+{Pi) 

In deriving Eq. Hll|) an explicit symmetry between the positive and negative charges has been assumed, i.e., iV+ — 

iV_|_+ = N , and = N |_. The first term in this expression is proportional to the unpolluted balance function, 

B. The second term can be simplified by assuming that the surplus charges are uncorrelated with the other charges 
and that they are also uncorrelated with themselves, aside from overall conservation of charge. 

N+{Pi)NssiPi,P2) = Ns{Pi)Ns+iPi,P2)^^, (12) 
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where Q represents the maximum integrated surplus charge. The balance function can then be expressed as: 



' 2{N+{P^) + Ns{Pi)) ' ' ' Q-1 2(7V+(Pi)+iV5(Pi))' 
B{P2\Pi) ^ + ^ l^(p^) ■ 

If the charges used to construct the balance function obey strict charge conservation, a perfect detector would satisfy 
the normalization conditions, 

P2 

J2NsiP,)^Q, 
P2 

P2 ^ 

After inserting Eq. (|14() into Eq. H13() . one can see that the normalization of the balance function is unchanged by 
the surplus charge. However, the shape is altered as the balance function is comprised of two components. The first 
term in Eq. H13(l describes the separation of balancing charges, while the second term is governed by the separation 
of two random balancing charges. The relative weights of the two terms is determined by the fraction of the charge 
which owes itself to a surplus in the initial state. Thus, the effect of the surplus charge is to dampen the contribution 
from the balancing charges and to average in a second contribution. 

For the STAR detector at RHIC, this second term is fairly small even for protons as the number of surplus protons 
is less than 10 per unit rapidity in central collisions Given that charge conservation constraints would suggest 
Q = 158, the effect of the extra charge is to first dampen the balance function by approximately 15%, and secondly 
to add in a second component whose width is characteristic of the acceptance, and whose magnitude is only one or 
two percent of the contribution from balancing charges. 

In order to eliminate the contribution from surplus charge and determine B from experiment, one can consider an 
object similar as to what is used to create the balance function numerator using mixed events. This object will be 
referred to as M{Pi,P2) and will be constructed from mixed events, where N^j^{Pi,P2) signifies that the charges a 
and b which satisfy the momentum constraints Pi and P2 are chosen from separate events. 

M{Pi,P2) =iv™_(Pi, P2) - iv;'V(Pi, P2) + iV!"+(Pi, P2) - iv"_(Pi, P2) 

-iv;"_(Pi, P2) - iv^+(Pi, P2) + iV!"+(Pi, P2) - iV!"_(Pi, P2) (15) 

- Nl'siPuP2) - Nr+iPi,P2) - iVJ?(Pi, P2) + iV™,(Pi, P2) + 7V,"l(Pi, P2). 

Since the counts for different events are independent, Nl'_ = N'l^ = N""^ = N"'_, = TVJ'l and TV™^ = iV™^. 
Thus, M becomes 

M(Pi,P2) = -iV^S(Pi,P2). (16) 

One could define a similar object using pairs from the same event, 

7V(Pi, P2) =Ar+_(Pi, P2) - 7V++(Pi, P2) + iV_+(Pi, P2) - iV— (Pi, P2) 

=7V+_(Pi, P2) - 7V++(Pi, P2) + iV_+(Pi, P2) - 7V__(Pi, P2) - iV55(Pi, P2). 

Again if one assumes that the only correlation between surplus charges involves a constraint on the overall number 
of surplus charges, 

iV(Pi,P2) = 7V+-(Pl,P2) -iV++(Pi,P2) +7V_+(Pi,P2) -7V__(Pl,P2) - ^^iVj?(Pi,P2). (18) 

One can then generate the desired balance function i?(P2|Pi) by combining N and AI, 

-, , , N(Pi,P2)- M(Pi,P2)(Q-l)/Q , , 

B(P2 Pi) = ^ ' ' \ \ '^^ (19) 

V 21 17 7V< Pi 
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where iV< is if the surplus charges are positive and is if the surplus charges are negative. 

It should be emphasized that this derivation assumed that the detector has equal acceptance for positive and 
negative charges. The ability of B to ignore the polluting surplus charge is based on the assumption that the surplus 
charges are uncorrelated with one another and are equally correlated with the "+" and "-" charges. As discussed in 
the next section and in 0, such correlations can be important, especially at small relative momentum. 

To illustrate the importance of these corrections, we generate a pp balance function from a simple model of a 
boost-invariant emission of particles governed by a temperature of 120 MeV and a maximum transverse velocity of 

0. 7c. It is assumed that the number of protons per unit rapidity is 28 and that the number of antiprotons is 21 to be 
consistent with measurements from RHIC The polluted balance function as described in Eq. H13|) is displayed 
in Fig. along with the corrected balance function B. This calculation is generated by assuming that particles were 
emitted from sources with random rapidities, but that two balancing particles are emitted from sources with the same 
velocity. The parameter Q used in Eq. ((T^ is assumed to be 158. 

Two additional modifications have been added to Eq. H13I) in order to more fairly illustrate the magnitude of the 
effect of the surplus charge. First, the function B was scaled down by 40% to account for the fact that the charge 
of an antiproton is often balanced by a neutron or by a A. Secondly, the simulated momenta were put through an 
acceptance filter which crudely mocks the acceptance of the STAR detector at RHIC. Particles were required to have a 
Pt greater than 100 MeV/c and a momentum of magnitude less than 700 MeV/c. The pseudo-rapidities were confined 
to a region of midrapidity, —1.1 < ry < 1.1. 

As illustrated in Fig. O the effects of the extra charge are mainly to dampen the balance function. The importance 
of correcting for the surplus charge would certainly be magnified if one were to analyze balance functions from SPS or 
AGS collisions where the fraction of extra protons is much higher. These corrections are not-model dependent, and 
the corrected balance functions exactly reproduce B. However, it should be emphasized that this statements relies 
on the assumption that the surplus charge is uncorrelated with other surplus charges, and with the pair- wise created 
charges. 

IV. FINAL-STATE INTERACTION DISTORTION TO THE BALANCE FUNCTION 

The balance function is implicitly predicated on the assumption that there are no residual correlations between a 
given charge and all other charges besides its balancing partner, i.e., all other charges are statistically eliminated from 
the distribution by the like-sign subtraction. Not all correlations cause problems. For instance, flow correlations tend 
to be identical between particles of the opposite charge or the same charge and thus fall out of the balance function. 
On the other hand final-state interactions involve all the other charges and depend sensitively on the relative signs 
of the charges. This distortion can rise linearly with the multiplicity since the number of charges with which a given 
charge can correlate rises linearly with the multiplicity. However, correlation functions tend to approach unity at 
higher multiplicity in accordance with expectations for increasing source size. This makes the resulting multiplicity 
dependence of the distortion non-trivial. 

A method for estimating the distortion to the balance function from residual interactions was provided in "^l. The 
same method is applied here. For every balancing pair, pa and pb, one must consider the correlation weight with other 
pairs whose momenta are Pc and pd- The weight, w{pa,Pb'iPc,Pd) can be estimated: 

w{Pa,Pb;Pc,Pd) ~ C++{pa,Pc)C {pb,Pd)C+-{pa,Pd)C-+{pb,Pc)- (20) 

Ideally, the balance function would isolate the ab pair and the interaction with the cd pair would cancel from the 
subtraction, — A^++. The correlations will lead to distortions if 

w{Pa,Pb;Pc,Pd) ^ w{pa,Pb;Pd,Pc)- (21) 

1. e., distortions are caused only by those interactions which differ between same-sign and opposite-sign particles. 
For instance, an isoscalar exchange of between pions would not bring on a distortion, but a Coulomb interaction or 
identical-particle interference would provide a source for distortion. 

We simulate these effects for tt+tt^ balance functions with the same blast wave model described in Section ^1 
The source of the ab pairs was chosen to move with one randomly chosen velocity, while the cd source was chosen 
to move with a different velocity. The distributions were calculated for the balance function numerators using the 
four particles, but rather than incrementing he distributions by unity, the distributions were incremented by the 
weight described in Eq. H20|) . By considering the contributions due to the extra weight separately from the usual 
contributions arising between the balancing particles, and weighting them appropriately for the given value of dn/dy, 
balance functions were calculated with and without the distortions. 
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In addition to the blast-wave parameters, the residual-interaction distortion is sensitive to the form for the two 
particle correlations and to the pion multiplicity. The multiplicity of charged pions was chosen to be 300 for both 
positive and negative pions, roughly consistent with measurements at RHIC |l8l |. The correlation functions were 
chosen to correspond to Gaussian sources of radius, i?i„v = 6.0 fm, again consistent with measurements at RHIC |19i]. 
A more sophisticated treatment would account for the pt dependence and the directional dependence of the radius. The 
correlation calculation presented here accounted for both identical-particle interference and for the mutual Coulomb 
interaction. 

The results of the calculation are shown in Fig.sQlandEl A crude filter for the STAR detector at RHIC was applied 
and it was assumed that 70% of the pions had their charge balanced by other pions, which affects the normalization 
of the balance function. 

Residual interactions can either strengthen or diminish the balance function depending on the relative momentum. 
At very small relative momentum, the balance function at small relative momentum rises due to the Coulomb en- 
hancement of the TT+TT^ correlation function and the Coulomb repulsion of the tt+tt-I- and 7r~7r~ correlation functions. 
For values of Qinv larger than a few MeV/c but less than ^ 25 MeV/c, the identical-particle interference which 
enhances the tt'^tt'^ and tt~tt~ correlation functions, diminishes the balance function since same-sign pairs contribute 
negatively to the balance function. At larger relative momenta, Coulomb effects again dominate. The effects are less 
dramatic when the balance function is viewed as a function of relative rapidity. 

The distortion of the balance function in Fig. ^ is dominated by Coulomb effects at large momentum. The 
correlation weights are driven by the squared quantum wave function. However, the correlation for large values of 
qR can be understood by considering the classical analog to the wave function. As shown in Rcf. the classical 
analog is: 



d^qi _ qi_ 

(22) 



4>^iq,r) 

u~ 

'l - 8ZqZban/rqj 
wl - AZaZba^i/rqj, 

where a is the fine structure constant, the product of the charges of the two species is ZaZh, fi is the reduced mass 
and qi and qf = Qinv are the initial and final relative momenta. Thus, the effects of Coulomb interactions, in the 
classical limit, only diminish as a function of XjQl^^. By averaging 1/r over a Gaussian source characterized by the 
Gaussian source size i?, one can find the asymptotic form for the classical correlation function. 

Cclass('3inv) ~ 1 7^9 (23) 

The classical result for the tt+tt" correlation function is compared to the quantum result for a six-fm source in Fig. 
El The agreement is remarkable for Qinv > 25 MeV/c when qR > 1, especially for the opposite-sign case where there 
is no identical-particle interference. 

Even though the correlation function goes to zero proportional to l/Qf^^, the phase space is increasing as Qf^^- 
Thus, the Coulomb interaction remains important to remarkably large momenta. For more central collisions, the 
value of R rises, but the number of particles with which a given particle is correlated also rises. If the multiplicity 
scales as R^, it is clear that the Coulomb distortion will become acute for central collisions. 

Strong-interaction distortions have not been considered in these calculations. In terms of the tttt phase shifts, Si, 
the contribution to the correlation function from strong interactions can be approximated by the relation |2lj |. 

The strength of the strong-interaction correlation falls much more quickly with R than does Coulomb-induced corre- 
lation. If R^ were to scale linearly with multiplicity, the effect of the strong interaction on the balance function would 
be roughly independent of multiplicity or centrality. This difference in the behavior derives from the fact that a given 
pion interacts with only its neighbors through the strong interaction while it may interact with nearly all particles 
through the Coulomb interaction. If the breakup density is independent of centrality, the number of neighbors stays 
constant and the distortion to the balance function from strong interactions should not be strongly centrality depen- 
dent. On the other hand, the Coulomb distortion interaction should be much stronger for central collisions than for 
peripheral collisions. 

The ingredients for calculating the distortion were the correlation weights and spectra along with the procedure for 
generating the pairs, {pa,Pb) and (pcPd)- In principal, the spectra and correlation weights can be taken or inferred 
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from data without introducing a theoretical modcL However, the generation of the pairs can not be extracted directly 
from data due to the correlations between pa and ph and those between pc and pd- Since particles are produced 
pairwise, it is necessary to include these correlations because the inter-pair interaction must attract pairs rather than 
single particles. I.e., the net charge in the medium can not change, but it can be polarized. At face value, this is 
an explicit model dependence. However, the parameters that govern the correlation between a and h and between 
c and d are precisely those parameters used to model the undistorted balance function. Thus, no additional model 
parameters would be introduced to calculate the distortion. Thus, the distortion from residual interactions can not 
be subtracted from experimental results in a model-independent fashion, but it can be modeled theoretically without 
additional parameters. 

Given the significant effects from inter-pair correlations, it is imperative that the balance function analyses correct for 
these distortions. Fortunately, the corrections can be confidently modeled, and the robustness of the balance function 
is not compromised. However, this conclusion is predicated on an understanding of the two-particle correlations. 
Since a correlation of a fraction of a percent can significantly alter the balance function, the issue of strong-interaction 
corrections to the balance functions should be revisited. 

Strong-interaction effects can be divided into two categories. The first category would be s-channel interaction 
which have particle-antiparticle channels, e.g. p'^ tt^tt". But, this source should not be considered as a distortion 
since the two pions are indeed a balancing pair. For instance, if all pions resulted from decays, the balance function 
would peak at the invariant mass of the p, and provide an important clue as to the creation mechanism for pions. Such 
resonant contributions can be calculated in a microscopic model or in a thermal calculation based on the canonical 
ensemble. A second source of strong-interaction effects is the interaction with other bodies through non-resonant 
interactions. Since the strong interaction is short range, this interaction should involve only a few neighbors. For 
large sources, the Coulomb interaction provides a larger effect on two-particle correlation function than does the 
strong interaction. Nonetheless, it would be worthwhile to better quantify the significance or insignificance of the 
strong interaction. 

V. SUMMARY AND DISCUSSION 

Charge balance functions were developed with the hope of identifying balancing charges on a statistical basis. Two 
effects prevent the like-sign subtraction from accomplishing this goal to high precision. As shown in Sec. IIIII the 
excess nucleons coming from the colliding nuclei provide a modest pollution to the balance function in measurements 
at mid-rapidity at RHIC. Fortunately, these effects can be easily subtracted. The second source of distortion derives 
from the inter-pair interaction of balancing charges with other particles in the medium. These distortions become more 
important in high multiplicity events. As shown in Sec. IIVI for high multiplicity events these distortions are most 
strongly affected by the Coulomb interaction. These effects are also more noticeable for balance functions calculated 
in Qinv than they are for balance functions calculated in relative rapidity. Although it is difficult to subtract these 
distortions in a model- independent fashion, it is straight-forward to include these effects in a theoretical treatment. 
In addition to the typical parameters one would use to model balance functions, modeling the distortion requires only 
an additional understanding of two-particle correlations. As these correlations can be extracted from measurement, 
the distortion from inter-pair interactions can be modeled quite confidently. In central collisions, this distortion can 
be a 20% effect, and if the correlation functions are understood to the 90% level, the residual systematic uncertainty 
is probably of the order of one or two percent. 
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FIG. 2: The width of the balance function is shown for the three rnornenturn components. The calculations assumed a blast 
wave scenario with the collective velocities of the source points for balancing pions being separated longitudinally according to 
a Gaussian distribution of width u^. The calculations assumed a breakup temperature of 120 MeV and a maximum transverse 
collective velocity of 0.7c. 
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FIG. 3: Proton-antiproton balance functions are shown for a blast wave model with and without corrections for the surplus 
unbalanced protons. The corrected balance function (circles) is constructed assuming a breakup temperature of 120 MeV, a 
maximum transverse velocity of 0.7c, and cr,, = 0. The calculation was scaled down by 40% to account for balancing of the 
proton's charges by other species. The distorted balance function (squares) is based on the proton excess as measured by the 
BRAHMS collaboration. Both balance functions were filtered through the STAR acceptance. 
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FIG. 4: As a function of Qinv, 7r'''7r~ balance functions from a blast-wave model are shown with (circles) and without (squares) 
the distorting effects of inter-pair interactions. The model assumed a breakup temperature of 120 MeV, a maximum transverse 
velocity of 0.7c, and cr^ — 0. The undistorted balance function was scaled by 70% to account for balancing by other species, and 
both balance functions were filtered by the STAR acceptance. The significant enhancement for momenta between 60 MeV/c 
and 400 MeV/c owes itself to the Coulomb interaction between pions. 




FIG. 5: The same as Fig. |1] only with the balance function being plotted as a function of relative rapidity. The distorting 
effects are less noticeable in Ay. 
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FIG. 6: Correlations for same-sign (squares) and opposite-sign (circles) pions are shown for a Gaussian source of size R — 6 
fm in the upper panel. To illustrate the Coulomb effects that force C{Q) to approach unity as , C{Q) — 1 is multiplied 
by and displayed in the lower panel. The lines represent the constant expected for Coulomb interactions described in Eq. 
1231 1. Since phase space increases as Q^, inter-pair correlations distort the balance functions for relative momenta of several 
hundred MeV/c. 



